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1. Introduction 


The phenomena of flow separation and subsequent reattach- 
ment which occur due to a sudden expansion in the flow passages 
such as backward facing steps have been recognized as impor- 
tant industrial situations. This complex flow structure present in 
heating or cooling applications such as cooling electronic equip- 
ments, cooling turbines blades, combustion chambers, chemical 
processes, cooling of nuclear reactors, wide angle diffusers, high 
performance heat exchangers, energy systems equipment, and flow 
in valves. In many instances separation of flow is undesirable and 
leads to unwanted pressure drops and energy losses which require 
additional fan or pumping power to overcome them. However, in 
other circumstances flow separation may be encouraged, such as 
in burner flame stabilization use for turbulence promotion leading 
to enhanced mixing or heat and mass transfer rates. 

The problem of flow over backward facing step has been 
extensively investigated, both experimentally and numerically. 
However, case studies such as injection flow, ribs geometry, double 
step ducts, and onset inlet flow, are not considered in the present 
review. The current review will only consider the studies of the 
case of steady-state convective flow and heat transfer over single 
backward facing steps. 

A large number of experimental and numerical studies focus 
on flow and heat transfer behavior of convective flow over back- 
ward facing step geometry have been reported. These vast results 
and information have dealt with different conditions, parameters, 
geometry dimensions, and instrumentation, which indeed to unde- 
fined solid base for comparison purposes to indicate more accurate 
methodology for solving the case studies. 

The dispersion of these results were attracted the attention of 
many researchers to unify the information to a general criteria. 
In 1992 and 1993, the Aerospace Heat Transfer Committee (K- 
12) of the Heat Transfer Division of the ASME held a technical 
session for benchmark heat transfer problems during the Winter 
Annual Meeting to solve numerically laminar mixed convection 
flow over two-dimensional horizontal and vertical backward fac- 
ing step, respectively [1-3]. Moreover, due to the importance of the 
separation and reattachment phenomenon Abu-Mulaweh in 2003 
[4] reviewed the results of flow and heat transfer of single-phase 
laminar mixed convection flow over different orientations of both 
backward and forward facing steps for several previous work. Most 
ofthe studies in this literature declared that only one large primary 
recirculation region develops downstream of the backward-facing 
step, while in some cases a secondary recirculation may develop 
between the downstream stepped wall and the step or at the top 
wall of the duct due to increase in step height and inlet velocity. 

The purpose of this paper is threefold. The primary purpose is to 
review the available results of the flow and heat transfer presented 
in previously published data for single backward facing step. The 
results of massive number of different parameters affect the fluid 
flow and heat transfer characteristics are summarized and pre- 
sented with a large number of correlation equations to predict the 
peaks of the flow and heat transfer in the recirculation regions. The 


Nomenclature 

AR aspect ratio (W/s) 

Cp specific heat (kJ/kg K) 

ER expansion ratio (H/h) 

f elliptic relaxation function 

Gr Grashof number (gB ATs? [12) 

g gravitational acceleration 

H total channel height (m) 

h inlet channel height (m) 

k thermal conductivity (W/m K) 

k turbulent kinetic energy (m?/s?) 

L total length of the channel (m) 

Nu local Nusselt number (hL/k) 

Nu max local Nusselt number (hL/k) 

P pressure (Pa) 

q heat flux (W/m?) 

Re Reynolds number (pUoh/ u) 

Ri Richardson number 

S step height (m) 

T temperature (K) 

Uo bulk velocity at the inlet (m/s) 

u velocity component x direction (m/s) 
v velocity component y direction (m/s) 
w velocity component z direction (m/s) 


Xuw-lines [(du/dy)z=constlly=o] 

X,Y,Z Coordinate directions 

Xi, Xe upstream, downstream lengths (m) 
Xr, Xo, Xn peak points 


Greek symbols 

thermal expansion coefficient (1/T) 
density (kg/m?) 

dynamic viscosity (N m/s) 
nanoparticles volume fraction 
amount of difference 

inclination angle 

kinematic viscosity (m2/s) 
dimensionless oscillation frequency 
turbulent dissipation rate (m?/s3) 
buoyancy parameter (Gr;[Re;?) 
v72[k 


wir 7g cS PSE 


Subscripts 
0 inlet 
w wall 


second purpose is to understand the characteristics and functions 
of nanofluids, to expect their effects and heat transfer enhance- 
ment in such geometries. The third purpose is to open a gate for 
new researches and propose suggestions that could lead to improve 
the ability to predict separation and reattachment phenomenon 
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using nanofluids. The review starts with an extensive review on 
the flow and heat transfer in backward facing step channels. After 
that, a comprehensive review for nanofluids and its characteristics 
is described. At the end, the review concentrates on the flow and 
heat transfer over a backward facing step using nanofluids. 


2. Flow geometry 


The general feature of backward facing step flow has the sim- 
plest reattachment flow, but the flow field is still very complex. 
The scope of this review is to summarize the previous efforts 
that studied the single-phase convective flow over two and three 
dimensional horizontal single backward facing step expansion in 
horizontal, inclined, and vertical orientations. A schematic diagram 
for backward facing step at either uniform wall temperature (UWT) 
or uniform heat flux (UHF) is shown in Fig. 1. 

The advent of new instrumentation such as laser anemometer, 
pulsed-wire anemometer, laser-Doppler measurement, hot wire 
probes, etc., have made a renaissance for proliferation of new 
researches in this area and facilitate the study of sophisticated 
three-dimensional flow structures associated with separation and 
reattachment phenomena in this geometry. In experimental stud- 
ies, the fluid velocity and temperature distributions were measured 
using Laser-Doppler velocitimeter and cold-wire anemometer, 
respectively, at any location. In numerical studies, the step geom- 
etry and its boundary conditions were modeled using Boussinesq 
approximation. The steady three-dimensional mass conservation, 
momentum equations, and energy equation governing the fluid 
flow and heat transfer problem are reduced to the following forms 
[5]: 


Continuity equation: 
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flow over a backward-facing step was done in the late 1950s. The 
advent of new instrumentation and the improvement of the numer- 
ical codes increase the number of new researches in such problem 
and facilitate the complex study of three-dimensional flow in the 
separation and reattachment zone. The horizontal, inclined, ver- 
tical cases were investigated for different geometrical, boundary 
conditions and fluids properties. 


3.1. Horizontal flow 


One of the first investigations dealt with laminar regime and 
subsonic free stream flowing over a step was presented by Gold- 
stein et al. [6]. They concluded, based on their experimental 
measurements, that in laminar regime the reattachment point is 
not a constant value as in turbulent regime, but it depends on the 
boundary displacement thickness and the hydraulic Reynolds num- 
ber. Interest in studying the reattachment phenomena over the 
backward-facing step continued in the 1970s and 1980s. Experi- 
mental and numerical studies for 2D pure forced convection flow 
where there is no heat transfer effect have been investigated by 
Denham and Patrick [7], and Armaly et al. [8,9]. 

Denham and Patrick [7] studied the measurements of 2D lam- 
inar flow over backward-facing step. The water flow in a range 
of Reynolds number (50 < Re < 250) passed through a single plane 
duct expansion with an expansion ratio and aspect ratio of ER=3, 
and AR ~ 20, respectively. They indicated that the general flow field 
bears similarities to other recirculating flows like axi-symmetric 
and the 2D duct sudden expansions but recirculation zone lengths 
and recirculated mass flow rates were smaller at a given value of 
Reynolds number. They noticed that a small fluctuation occurs at 
Re=229 while the other lower values remain steady; this fluctu- 
ation observed to be the onset of transition to turbulence of the 
separated boundary layers. Furthermore, it was found that there 
is a tendency towards secondary separation at the upper wall as 
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Z-momentum equation: 
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Energy equation: 
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There are different techniques were used to solve the partial dif- 
ferential equations along with specific boundary conditions, such 
as finite difference (FD), finite volume (FV), finite element (FE), or 
lattice Boltzmann (LB). These methods are solved by using differ- 
ent numerical software codes such as FLOTRAN, FLUENT, TEACH, 
COUGAR, KAMELEONII, FIDAP, NEKTON, CUTEFLOWS, QMRCGS, 
etc. 


3. Backward facing step geometry 


Convective flow and heat transfer over a backward facing step 
has been extensively investigated in the last decades, both exper- 
imentally and numerically. It has been reported in the literature 
that the first effort for studying the separation and reattachment 


reattachment occurs, this phenomenon was produced due to the 
inflections of the velocity profiles and it was clearer at the highest 
Reynolds number. 

The effect of Reynolds number for different flow regimes was 
conducted by Armaly et al. [8]. They reported the measure- 
ments of velocity distribution and reattachment length through 
downstream of a single backward-facing step mounted in a two- 
dimensional channel. Air was used as a working medium with 
Reynolds number range of 70 « Re « 8000 in a test section has a 
high aspect ratio of AR- 36. Numerical predictions were also pre- 
sented to be compared with the experimental results as the flow 
maintained its two-dimensionality in the experiments. The com- 
parison showed a good agreement with the experimental findings 
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Fig. 1. Schematic diagrams for backward-facing step. 


up to a Reynolds number Re« 400. Armaly et al. [9] showed that 
the vortices start to develop in the laminar region at Re~ 800, 
but it becomes very pronounced at Re - 2300, and damped out at 
Re=5000. It was indicated that in both the pure laminar (Re < 400) 
and the fully turbulent (Re = 6600) flow regimes, the flow is indeed 
two-dimensional over large parts of the width of the test section. 

Most of the existing research efforts have been focused on two- 
dimensional situations, though most of the realistic applications 
are practically dominated by three dimensional geometries. In the 
early 1990s many researchers tend to study the flow field in three- 
dimensional due to the presence of certain effects which are not 
visualized in two-dimensional (at centerline of the plane) such as 
the side-wall effect. The three-dimensional pure forced convection 
flow have been extensively investigated by Shih and Ho [10], Chi- 
ang et al. [11], De Brederode and Bradshaw [12], Hertzberg and Ho 
[13], Tylli et al. [14], Armaly et al. [15], and Nie and Armaly [16]. 
Shih and Ho [10] examined the water flow field behind a small 
aspect ratio backward-facing step with AR -3 using Laser-Doppler 
Anemometer. The measurements ofthe three velocity components 
inside the separation region indicated that the flow is highly three- 
dimensional. The mean reattachment line is highly distorted across 
thespan with the shortest reattachment position located at the cen- 
treline and about 20% shorter than the value closer to the side wall. 
The maximum spanwise velocity component was indicated to be 
of the same order as the streamwise velocity at the point near the 
center of the step. 

Chiang et al. [11] aimed to clarify and present a deep expla- 
nation on skin-friction fields as a guide for exploring the physical 
interpretation of the three-dimensional flow separation behind the 
step, reattachment of the flow, and the roof recirculatory flow pat- 
tern. The channel designed to have an expansion ratio ER = 1.9432 
where the step was 0.9432. The parametric study was carried out 
by varying Reynolds numbers and various spans with values not 


greater than 800 and 10, respectively. It was shown that the flow 
separation-attachment occurs on the roof of the channel but it con- 
fined only to the sidewall region. The roof recirculation zone size is 
dependent on Reynolds number. 

The distance between the side-walls plays a significant effect 
with aspect ratio of the backward-facing step ducts. De Brederode 
and Bradshaw [12] systematically studied the side wall effects 
with top hat free stream velocity profile. It was concluded that the 
three-dimensional features could be neglected if the aspect ratio of 
the step was larger than 10. Furthermore, Hertzberg and Ho [13] 
observed that the presence of the strong spanwise velocity in a 
small aspect ratio rectangular sudden expansion will obviate the 
formation of the recirculation zone after the step and therefore has 
a zero reattachment length. 

Tylli et al. [14] studied the sidewall effects in three-dimensional 
laminar water flow over a backward-facing step experimen- 
tally and numerically. Experiments were based on both digital 
and stereoscopic particle image velocimetry; a spectral element 
method was used for the laminar, transition, and turbulent simula- 
tions. The geometry has an expansion ratio, ER, and aspect ratio, AR, 
of 2 and 20, respectively. It was noticed that excellent agreement 
between the 3D simulation and the experiments is obtained when 
Reynolds number reaches Rez 650. Moreover, it was found that 
at low Reynolds numbers (Re « 400), side wall effects do not affect 
the structure of laminar flow in the channel mid plane; thus, the 
mid plane field is accurately predicted by two-dimensional simu- 
lations. In addition, it was found that at higher Reynolds numbers, 
laminar flow is characterized by sidewall separation and the for- 
mation of a recirculation zone, which, however, does not penetrate 
up to the channel mid plane. The wall-jet intensity and flow three- 
dimensionality increase with Reynolds number for laminar flow, 
and decrease for transitional flow; three-dimensional variation is 
hardly noticed in the time average field of turbulent flow. 
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Armaly et al. [15] measured the velocity for three-dimensional 
laminar separated airflow adjacent to a backward-facing step using 
two-component laser Doppler velocimeter. The backward-facing 
step, with a height of s- 1.0 cm has an aspect ratio of AR=8, and 
an expansion ratio of ER - 2.02. The flow measurements were cov- 
ered Reynolds number range between 98.5 < Re < 525 to develop 
a clear view of both sidewall and the step effects, and compared 
with the previous 2D studies. It was determined that for Re < 98.5 
there was no recirculation flow region adjacent to the sidewall and 
at Re < 190, a small recirculation flow region was detected in the 
upper corner of the sidewall. It was concluded that the spanwise 
and the transverse peak velocities move to the center of the duct 
as the distance from the step increases. 

Nie and Armaly [16] measured the velocity adjacent to the 
bounding walls of three-dimensional backward-facing step flow 
using Laser-Doppler velocimeter. The study was performed for the 
purpose of mapping the boundaries of the reverse flow regions that 
develop adjacent to the walls of a geometry with a step height 
s=1cm, aspect ratio AR=8, and an expansion ratio ER- 2.02 as a 
function of Reynolds number. The Reynolds numbers were pre- 
sented in the range of 100 « Re « 8000 to cover the three types of 
flow regimes. The boundaries ofthe reverse flow regions were iden- 
tified by locating the streamwise coordinates on a plane adjacent 
to the bounding walls where the mean streamwise velocity com- 
ponent is zero. The three-dimensional reattachment point at the 
stepped wall was compared with the two-dimensional at the cen- 
ter ofthe test section. It was noticed that the three-dimension flow 
results are slightly higher in the laminar flow regime (Re « 400); 
significantly lower in the transition flow regime (400 « Re « 3400); 
and slightly lower in the fully turbulent flow regime as (Re > 3400). 
It was noticed that the size ofthe reverse flow regions increases and 
moves further downstream in the laminar flow regime; decreases 
and moves upstream in the transitional flow regime; and remains 
almost constant or diminishes in the turbulent flow regime; as the 
Reynolds number increases. 

The literature review reveals that there were several studies 
focused on the heat transfer phenomena over the backward-facing 
step which exhibit an influence in the fluids flow behavior. Aung 
[17], and Sparrow and Chuck [18] are the first investigators who 
established and reported an experimental and numerical results 
on heat transfer for 2D laminar air flow passing a backward-facing 
step channel heating at uniform temperature from below. In their 
studies, they outlined that for laminar airflow the heat transfer 
increases monotonically in the stream-wise direction and quan- 
titatively it is less than that of the flat plate value. It was concluded 
that the Nusselt number was independent of the Reynolds number; 
the local Nusselt number distribution begins with a low value at 
the step, and then increases monotonically and attains a maximum 
value at a position near the reattachment point then it decreases 
monotonically towards the fully developed value. 

Khanafer et al. [19] analyzed mixed convection of laminar pul- 
satile flow and heat transfer past a backward-facing step in a 
channel. Fluid flow and heat transfer characteristics were exam- 
ined in the domain of the Reynolds number, Richardson number 
and the dimensionless oscillation frequency as: 100 < Re < 1000, 
1.78 x 10-3 «Ri «10, and 0.1<@ «5. It was found that the 
elimination of the separation region occurs for Ri»0.1, and as 
Reynolds number increases, the recirculation zone along the 
heated downstream surface becomes larger. The thermal bound- 
ary layer decreases in thickness as Reynolds number increases and 
consequently increases the heat transfer rate. The average Nus- 
selt number and the length of the recirculation zone decrease 
with an increase in Richardson number. It was also noticed 
that the dimensionless local variation of skin friction coefficient 
along the heated wall increases with increasing the Richardson 
number. 


Chen et al. [20] simulated numerically the turbulent forced 
convective flow adjacent to a two-dimensional backward- 
facing step. The objective of this study was to explore 
the effects of step height on turbulent separated flow and 
heat transfer. The stepped wall downstream from the step 
was set to be uniform and constant heat flux qw -270W/m?, 
while other walls were treated as adiabatic. Reynolds number 
and duct's height downstream from the step were kept constant 
at Reg - 28,000 and H - 0.19 m, respectively. On the other hand, 
the expansion ratio was maintained to ER» 1.11, 1.25, and 1.67, 
respectively. It was noticed that as the step height increases; the 
primary and secondary recirculation regions, magnitude of the 
maximum turbulent kinetic energy increase, and the bulk temper- 
ature increases rapidly as well. Near the step and below the step 
height, it was found that the kinetic energy becomes smaller as the 
step height increases. Moreover, it was found that the maximum 
temperature becomes greater as the step height increases. The peak 
values of the transverse velocity component become smaller as the 
step height increases, but inside the recirculation region; the skin 
friction coefficient becomes less significant and smaller in magni- 
tude with the increase of the step height. 

The effects of the flow and heat transfer in a 3D backward-facing 
step ducts have been extensively studied using different parame- 
ters by Iwai et al. [21], Nie and Armaly [22], Armaly et al. [23], 
Barbosa Saldana and Anand [24], Barbosa Saldana et al. [5], and 
Lan et al. [25]. Iwai et al. [21] performed a numerical simulation 
for three-dimensional mixed convection flows over a backward- 
facing step at low Reynolds number in order to investigate the 
effects of the duct aspect ratio. Numerical results were conducted 
for an aspect ratio range between 4 < AR < 24, and Reynolds num- 
ber range between 125 < Re < 375 to study its effects on Nusselt 
number and the skin friction coefficient on the bottom wall. It was 
found that an aspect ratio of 16 at least is needed to secure 2D 
region in the central part of the duct at a Reynolds number of 250. 
It was also found that the maximum Nusselt number increases with 
an increase of the aspect ratio and Reynolds number. The skin fric- 
tion coefficient is found to be an analogous to the case of Nusselt 
number distribution. 

Nie and Armaly [22] presented a numerical analysis for a 
three-dimensional incompressible laminar forced convection flow 
adjacent to backward-facing step in rectangular duct. The effects 
of step height on the flow and heat transfer characteristics were 
the main objectives of this study. The geometry, thermophysical 
properties, and flow condition are based on the available measure- 
ments of Armaly et al. [23]. The simulations were performed at 
Re- 343 for different step heights (s- 0.008, 0.010, and 0.012 m), 
while the other geometries were kept constant. It was found that 
increasing the step height increases the reattachment length, the 
Nusselt number, the size of the sidewall reverse flow region, and 
the general three-dimensional features of the flow. In addition, 
as the step height increases, the distance between the step and 
the point where the impingement flow reaches the stepped wall 
downstream increases, and a secondary recirculation flow region 
develops adjacent to the bottom corner of the step. For the case 
of small step height (s = 0.008 m), the minimum Nusselt number is 
located near the bottom corner between the step and the center- 
line of the duct. As the step height increases, the minimum Nusselt 
number is appeared to move towards near the bottom corner of 
the step and the sidewall. It was found that the friction coefficient 
increases in the streamwise direction along the centerline of the 
duct inside the primary recirculation flow region, but it decreases 
outside the recirculation region with increasing step height. 

Barbosa Saldana and Anand [24] studied numerically forced con- 
vective airflow over a three-dimensional backward-facing step. The 
horizontal backward-facing step channel was set to have a step 
height s» 1 cm, an aspect ratio of 8 and an expansion ratio of 2. The 
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stepped wall was fixed to a constant heat flux qw = 50 W/m? while 
the other walls were considered as insulated. The Reynolds num- 
bers were set to be in the range of 98.5 and 512. The local Nusselt 
number along the bottom wall was found to lie in the vicinity of the 
Xy-line and the xy-line points of intersection where the shear stress 
along the bottom wall is equal to zero. Furthermore, the velocity 
profiles revealed that for Reynolds number greater than 343 the 
flow does not reach fully developed conditions at the channel exit. 

Barbosa Saldana et al. [5] simulated numerically laminar mixed 
convective flow over a three-dimensional horizontal backward- 
facing step. All the walls are treated adiabatically except the bottom 
wall of the channel is subjected to a constant high temperature. 
An airflow of fixed Reynolds number Re - 200 enters hydrodynam- 
ically fully developed and isothermally in a duct with an aspect 
ratio of 4, expansion ratio of 2, and total length of 52 times the 
step height L= 52S. It was reported that as Richardson number (Ri) 
increases; the location of xy-line displaces to the upstream (higher 
velocity components), the recirculation zone is shortened, verti- 
cal size of the recirculation zone reduced, and negative values of u 
velocity (streamwise direction) next to the bottom wall ofthe chan- 
nel increases. In addition, as the Richardson number increased; the 
location of the maximum averaged Nusselt number distribution 
was moved and the maximum averaged streamwise shear stress 
distribution moved its location further upstream. 

Lan et al. [25] simulated the turbulent forced convection adja- 
cent to three-dimensional backward-facing step in a rectangular 
duct using a (k— € — £ — f) turbulent model. The aim of the study 
was to present a clear, realistic, and more accurate investigation 
near the walls to predict the turbulent heat transfer separated 
and wall-bounded flows better than the commonly used two- 
equation turbulence model. A User Defined Function (UDF) was 
implemented using FLUENT-CFD, and two-dimensional benchmark 
problems were simulated to validate the three-dimensional simu- 
lation code. The geometry has an expansion ratio and a step height 
of 1.48 and 4.8 mm, respectively. The three aspect ratios of 3, 8 and 
infinity (2-D simulation) were considered for studying its effect on 
the flow and heat transfer. The effect of Reynolds number in the 
range of 20,000 < Re < 50,000 was also considered. The effects of 
the Reynolds number and aspect ratio on the flow reattachment 
was minimal in the range of parameters that were examined, but 
they were shown to have significant influence on the heat trans- 
fer. The reattachment length was found to be constant in central 
portion of duct width, but starts to increase as it approaches the 
side walls and reaches its maximum value on the side walls. It was 
noticed in the primary recirculation flow region and in the vicin- 
ity of the recirculation line a significant spanwise flow developed. 
The maximum local Nusselt number was found to be developed in 
the region where the spanwise velocity is maximum adjacent to 
the heated wall close to the sidewall. The bulk air temperature was 
found to decrease as the Reynolds number increases, but it was not 
influenced significantly by the changes in the aspect ratio. 


3.2. Inclined flow 


The backward-facing step channel orientation effects attracted 
some researchers such as Lin et al. [26], Lin et al. [27], Hong et al. 
[28], Abu-Mulaweh et al. [29], and Iwai et al. [30] to investigate 
the separation region and the reattachment point with respect 
to other parameters. The literature declared that the first effort 
dealt with backward-facing step in inclined duct was introduced 
by Lin et al. [26]. They examined numerically the effects of inclina- 
tion angle (0° < $ < 90° from the vertical) on the two-dimensional 
mixed convective laminar flow and heat transfer in a duct with 
a backward-facing step under buoyancy-assisting flow conditions. 
The case study was taken from a previous study done by the same 


author for a vertical duct with a backward-facing step Lin et al. 
[27]. The geometry has an expansion ratio and step height, of 2 and 
0.48 cm, respectively. The straight wall at the top of the duct was 
maintained at the inlet air temperature of 20°C, the downstream 
stepped wall was maintained to a uniform temperature of 35°C, 
while the upstream wall and the step were maintained adiabatic. 
The Reynolds number was maintained at 50. It was found that the 
recirculation size increases as the orientation transfer from the ver- 
tical axis. The streamwise component of the Grashof number has 
a strong influence on the flow, but the transverse component has 
a small effect on the velocity and temperature profiles in the fully 
developed regime. The substantial fact elaborates that the reattach- 
ment length, fluid temperature, and the shear stress at the unheated 
wall increase as the inclination angle increases from the vertical. In 
contrary, for Nusselt number for both heated and unheated walls, 
the streamwise buoyancy force, and the shear stress at the heated 
wall are decreased. 

Hong et al. [28] were the first who investigated the buoyancy 
opposing effects in an inclined channel with a backward-facing 
step. They studied numerically the influence of the relevant inclina- 
tion angle and Prandtl number for two-dimensional laminar mixed 
convective airflow for buoyancy assisting and opposing flow condi- 
tions. The geometry has an expansion ratio and step height, of 2 and 
0.0048 m, respectively. The straight wall at the top of the duct was 
maintained at the inlet air temperature of 293 K, the downstream 
stepped wall was maintained to a uniform heat flux of 200 W/m2, 
while the upstream wall and the step were maintained adiabatic. 
The Reynolds number was fixed at 100, and the Grashof number 
was fixed at 609. The inclination angle was varied from 0° to 360°, 
while other parameters were kept constant. It was found that as the 
inclination angle increases from 0° to 180°, the reattachment length 
increases, but the wall friction coefficient and the Nusselt number 
at the heated wall decrease. Moreover, the transverse buoyancy 
force was noticed to have a less significant effect on the flow and 
heat transfer downstream the reattachment point, while it has a 
strong effect in the recirculation region. The Prandtl number was 
varied in the range of 0.07 < Pr < 100, while the inclination angle 
was fixed at 0°. It was found that as the Prandtl number increases, 
the Nusselt number and the reattachment length increase, but the 
wall friction coefficient decreases. 

The understanding of the separated zone and reattachment pro- 
cess was poor for the inclined channel with a backward-facing 
step due to the inability to predict simple reattaching flows over 
a wide range of parameters. The first experimental effort dealt 
with backward-facing step mounted in inclined channels is pre- 
sented by Abu-Mulaweh et al. [29]. They measured experimentally 
and simulated numerically the influence of buoyancy-assisting 
laminar mixed convection in boundary-layer flow over horizontal 
and inclined two-dimensional backward-facing step. The objec- 
tives of the study were to determine the reattachment length 
and the onset of vortex instability for different wall temperatures 
(0 x AT x 30°C), free stream velocities (0.285 <u. < 0.7 m/s), step 
heights (0.35 < s < 0.8 cm), and inclination angles (30° < @ < 90°) as 
measured from the vertical direction. The geometry of the test 
section was consisted of a heated constant temperature down- 
stream section, adiabatic upstream section and step, and the top 
section was a free stream. It was noticed in horizontal case that 
the buoyancy force has a negligible effect on the velocity and tem- 
perature distributions, but it significantly influences the onset of 
instability. On the other hand, it was found that as the inclina- 
tion angle increases from the vertical, the streamwise component 
of the buoyancy force decreases, which causes an increase in 
the velocity gradient and a decrease in the temperature gradient 
at the wall in the recirculation region. It was concluded that as 
the inclination angle increases from the vertical, the local Nus- 
selt number decreases, while the reattachment length and the 
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location of the maximum Nusselt number behind the backward- 
facing step increase. It was noticed that the maximum Nusselt 
number occurs downstream of the reattachment point. Moreover, 
the space between the maximum Nusselt number location and the 
reattachment decreases slowly as the inclination angle increases. 

As it has been introduced in the previous sections, all ofthe exist- 
ing researches for the inclined channel with a backward-facing step 
are focused on two-dimensional situations, though most of the real- 
istic applications are practically dominated by three dimensional 
geometries. It is declared from the literature that the first three- 
dimensional study is presented by Iwai et al. [30]. They simulated 
numerically the effect of inclination angle in three-dimensional 
mixed convective flow over a backward-facing step in a rectan- 
gular duct. The downstream wall was kept with a constant heat 
flux, while other walls were kept adiabatic. Parameters such as 
Reynolds number, expansion ratio, and aspect ratio were kept con- 
stant, at 125, 2, and 16, respectively. The effects of two inclination 
angles; $4 pitch angle changing in the range of 0° < $4 < 360°, and 
$» rolling angle changing in the range of 0? < $5 < 180° are the 
main interest of this study. The three positions ofthe reattachment 
point X;, the peak Nusselt number point Xy and the downstream 
end of the secondary recirculation region Xs on the centerline of 
the heated wall are shifted with $4. It was found that when 4 
is changed while keeping @2 to be zero degree, positions of the 
maximum Nusselt numbers are symmetrically obtained at the posi- 
tions on the heated wall near the two side walls, similar to the 
cases of pure forced convection. The maximum Nusselt number 
was appeared at the most upstream position in the case of $1 =0°, 
and the minimum at the most downstream position in the case of 
$41 = 180°. On the other hand, it was noticed that when $» is changed 
while keeping ¢; to be 90°, the flow and thermal fields become 
asymmetric about the duct centerline. It was observed that $» 
has little effect on the spanwise position of the maximum Nusselt 
number. 


3.3. Vertical flow 


It is indicated from the literature that the first effort of study- 
ing the effects of vertical backward-facing step was introduced 
by Lin et al. [31]. They examined and reported the numerical 
results of laminar mixed convection airflow and heat transfer, 
buoyancy-assisting through a two-dimensional vertical duct with 
a backward-facing step. The study was conducted to cover the 
domain from pure forced convective flow to the inlet starved 
convective flow where the average inlet velocity is smaller than 
the corresponding natural convective value. The geometry dimen- 
sions and boundary conditions were based on previous work done 
by Armaly et al. [7]. The influence of the temperature difference 
between the heated wall downstream of the step and free stream 
in the range of 1 < AT < 75°C was detected. It was found that as the 
wall temperature increases, the reattachment length decreases and 
the size of the secondary, inner recirculation region that develops at 
the corner between the heated wall and the step increases. On the 
other hand, the buoyancy level decreases for a fixed wall tempera- 
ture as the Reynolds number increases. It was noticed that the force 
of the pressure drop Fp is lower than the force of the kinetic energy 
F,, where a reverse flow could be developed from the downstream 
section of the duct along the cooler wall as outlined by Aung and 
Worku [32]. However, it was found that the disappearance of the 
reattachment from the heated wall will occur only in the starved 
flow regime. The Nusselt number was found to be increased as the 
wall temperature increases and the Reynolds number decreases. 
Both the peak Nusselt number length Xn, and the reattachment 
length of the main recirculation region, X;, decrease, while the 
length of the recirculation region, Xo, increases as the value of 
the buoyancy parameter increases. Moreover, the peak Nusselt 


number occurs downstream of the reattachment point and the 
spacing increases as the buoyancy parameter increases. 

The effects of step height in inclined ducts that presented before 
in this literature did not include the vertical orientation ($ - 0?). 
This gap of research motivated Hong et al. [4,33] to investigate 
the effects of other parameters due to different step height. They 
numerically studied the effects of the mixed convection flow and 
heat transfer through a two-dimensional vertical duct having a 
backward-facing step. The downstream wall from the step was 
subjected to a uniform heat flux (UHF), while the upstream wall 
and the step by itself were adiabatic, and the straight top wall was 
maintained at a uniform temperature that is equal to the inlet fluid 
temperature. Furthermore, they examined the effects of different 
ranges of parameters such as Grashof number, 0 < Gr; < 24,000, 
expansion ratio, 1.25 < ER < 4, and Reynolds number, 0 < Re < 150. 
It was found that Grashof number has a significant influence on 
the flow and heat transfer; when Grashof number increases, the 
wall friction coefficient and the local Nusselt number on the down- 
stream heated wall increase, and verse visa for the reattachment 
length. It was noticed that the reattachment length increases as 
the expansion ratio increases for less than 2.25, but for values 
greater than 2.25 the reattachment length decreases. Moreover, it 
was found that increasing the expansion ratio increases the fully 
developed value of the wall friction coefficient and the distance for 
the flow to reach the fully developed regime, but it decreases the 
peak and the fully developed Nusselt numbers. They have devel- 
oped correlation equations for reattachment length and the peak 
Nusselt number as shown in Table 1. 

Baek et al. [34] studied experimentally and numerically the 
buoyancy-assisting laminar mixed convection airflow and heat 
transfer over a vertical, two-dimensional backward-facing step. 
Different free-stream velocities, wall temperature differences, and 
step heights in the range of (0.37 < uo < 0.72 m/s), (10 < AT < 30°C), 
and (0.38 <s<1cm), respectively, were considered to visualize 
the flow and heat transfer behavior. It was found that for a fixed 
inlet velocity, the reattachment length, X;, decreases as the wall 
temperature increases. Similarly, for a fixed wall temperature, the 
reattachment length, X;, and the volume ofthe recirculation region 
decrease as the inlet velocity decreases. It was noticed that the 
recirculation region disappears and detaches from the heated wall 
and attaches only to the adiabatic step. Moreover, for the heat 
transfer visualization, it was found that the heat transfer rate is 
seen to increase as the temperature difference increases. The fluid 
temperature in the recirculation region behind the step (X«X;) 
was considerably higher than the temperature downstream of the 
reattachment point (X» X;). As the buoyancy force increases, the 
magnitude of the Nusselt number increases and the location of its 
peak moves closer to the step. An effort was made by the authors 
to modify the existing empirical correlation for the reattachment 
length in the laminar forced convection flow, as given by Goldstein 
et al. [6], to reflect also the influence of step height and buoyancy 
parameter. It was found that the location of the maximum Nus- 
selt number, Xp, and the reattachment length, X;, decreases as the 
buoyancy force increases. 

Abu-Mulaweh et al. [4,35] examined the influence of buoyancy- 
assisting laminar flow and heat transfer over a 2D vertical 
backward-facing step. The boundary conditions and general step 
geometry were taken from Baek et al. [34] except the downstream 
heated wall was maintained to a uniform heat flux in the range 
between 102 < qw « 290 W[m?. As the buoyancy force increases, 
the reattachment length Xr decreases but the local Nusselt num- 
ber increases and the location of its maximum values moves closer 
to the step. Moreover, it was found that the local Nusselt number to 
be lower in the region where X « X; and higher in the region where 
X» X; for higher free stream velocity. Correlation equations for the 
reattachment length were reported in Table 1. 


Table 1 
Summary of the research done on backward facing step. 
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Denham and 
Patrick [7] 


Armaly et al. [8] 


Shih and Ho [10] 


Chiang et al. [11] 

Tylli et al. [14] 

Armaly et al. [15] 

Nie and Armaly 
[16] 

Khanafer et al. 


[19] 


Chen et al. [20] 


Abu-Nada [169] 


Iwai et al. [21] 


Nie and Armaly 
[22] 

Saldana and 
Anand [24] 


Saldana et al. [5] 


Lan et al. [25] 


Lin et al. [26] 


Hong et al. [28] 
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NG 


50 < Re < 250 
ER=3 

AR«20 

70 « Re « 8000 
ER = 1.9423 
AR=36 

AR=3 

Laminar 
Re < 800 

ER = 1.9432 
AR=10.6022 
Lam < Re < Tur 
ER=2 

AR=20 

98.5 < Re < 525 
ER=2.02 
AR=8 

100 < Re < 8000 
ER=2.02 
AR=8 

100 < Re < 1000 
1.78E-3 <Ri<10 
0.1«co «0.5 
AR=8 

Re- 28000 
1.11 < ER x 1.67 
200 < Re < 600 
0<w «0.2 
ER-2 

125 < Re x 375 
4 

<AR<24 

Re =343 
8<s<12 mm 
98.5 < Re < 512 
ER=2 

AR=8 

Re- 200 

ER-2 

AR=4 

2E^ < Re < 5Ẹ4 
ER- 1.48 

3 <AR< œ 

0° < <90 
Re=50 

ER=2 

0° < ġ < 360° 
0.07 < Pr < 100 
Re=100 

Gr, =609 
ER=2 


X, =[5.462 — 1.62 cos(¢)] 
x [0.95 + 0.05 cos(29)][1 
+ 0.02 sin($)] 

X, 74.9 -1.7 
exp(—Pr/0.4) 


Xo = MAX(0.4462 cos($)|1 


— 0.5 sin(@)],0} 


Xo = 0.336 — 0.00165Pr 


Xn = [5.472 — 0.6535 cos($)][098 
* 0.02 cos(2¢)][1 + 0.02 sin()] 


Xn 7 3.53 * 1.71Pr, 
0.07 < Pr x 0.712 


Xn 74.753 — 0.0091Pr?, 


0.712 x Pr x 100 


Nus, max=[1.01 — 0.01 cos(29)] 
[1.752 * 0.151 cos($)] 
Nus,max = [0.52 + 1.29./Pr 
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Abu-Mulaweh 
et al. [29] 


Iwai et al. [30] 


Lin et al. [31] 


Hong et al. [33] 


Baek et al. [34] 


Abu-Mulaweh 
et al. [35] 


Abu-Mulaweh 
et al. [36] 


Iwai et al. [38] 


Abu-Mulaweh 
et al. [39] 


Abu-Mulaweh 
et al. [40] 


Inc 


Inc 


Ver? 


Ver 


Ver 


Ver 
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Ver 


UWT 


UHF 


UHF 
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UHF 
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UHF 
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Air 
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Air 
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Air 
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2D 


3D 


2D 


2D 


2D 


2D 
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3D 


2D 
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30° < ġ < 90° 

50 < Re, < 300 
3.5 <s < 8 mm 
0<AT<30°C 


0° <@, < 360° 
0° x4; < 180° 
Re=125 
ER=2, AR= 16 Ri=0.03 
$-0* 
Re, - 50, 100 
1<AT<75°C 
ER=2 

$-0* 

0 «Re, < 150 
0«Gr, < 2400 
1.25 «ER «4 


ó-0 
0.37 x us, x 0.72 
10x AT<30°C 
3.8<s<10mm 
p=0° 

50 < Re, < 330 
0<&<14E? 
3.5<s<8mm 
P= 180° 

0.41 <u < 0.66 
0<AT<30°C 
4<s<8mm 
$-0 
0<&<0.12 
Re=125 

ER=2, AR=16 
$-0 

O<u, x0.41 
S=22mm 

$-0 

Us; 70.41 m/s 
0<s<22 mm 


X, =0.023(8,/s) 

Re; + f(s[xi) - 0.023 
[5(xi/s)Re,;!? ]Res + f(s]x;) 
f(s[xi) 8995.8 (s|xj? — 114.2 
(s/xi) - 3.4 


Xr 7 Xo torcea (1 *(7.188 — 4.047R 
—1.758Re; EY 
Xrsforced = ao(R)* a1(R)Res 


ao(R) = — 1.252 + 1.345R — 0.1681R? 


a;(R) - — 0.02751 
0.04572R — 0.006414R? 
X, - (2.24 + 0.022Re,) 
exp[-42 &(s/x;) 986] 


Xr 7 Xrforced exp(—7.5 £2) 

X; force = 0.023(8;[s)Re; + f(s[xi) 
f(slxi) =8995.8 f(s|xi — 114.2 
(s/xi) — 3.4 


Nuüsmax = 1.658 + 0.0001483Gr," 
Nus max 74.862 — 2.172R 
+ 0.3292R? 


a Hor, horizontal; Inc, Inclined; Ver, vertical; UHF, uniform heat flux; UWT, uniform wall temperature; NF, nanofluids; E, experimental; N, numerical; #, not identified; -, no correlation equation. 
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As introduced in the previous section in this literature, the 
researchers were concentrated on buoyancy-assisting vertical 
flow. The first effort to study the influence of the buoyancy- 
opposing flow was conducted by Abu-Mulaweh et al. [36]. They 
performed experimental and numerical studies of laminar flow 
over a vertical backward-facing step to measure the effect of the 
buoyancy-opposing on the airflow and heat transfer behavior. The 
experimental setup, test section and backward-facing step geom- 
etry, and the boundary conditions, are identical to that used by 
Abu-Mulaweh et al. [35] except that the orientation of the air tunnel 
is changed by 180° to create buoyancy-opposing flow conditions. 
It was noticed for low buoyancy levels Ri « 4.4 x 1073; the length of 
the recirculation region downstream of the backward-facing step 
increases rapidly as the buoyancy level increases. However, for 
low buoyancy levels above 4.4 x 1073; the flow inside the recir- 
culation region remains laminar but the length of the recirculation 
region decreases rapidly due to the effect of the transition of the 
flow from laminar to turbulent and the flow becomes turbulent 
downstream of the recirculation region. The effect of temperature 
difference between the heated wall downstream the step and the 
free stream in the range of 0 < AT < 30°C was studied. It was found 
that as the temperature difference AT increases, the reattachment 
length increases and the Nusselt number decreases. 

Abu-Mulaweh et al. [4,37] studied the effects of heating the 
upstream wall and the step of two-dimensional, vertical backward- 
facing step geometry on the flow and heat transfer characteristics 
in the region downstream of the step. Four different upstream wall 
and step heating conditions were examined. The first case Tem- 
perature Boundary Condition (TBC) was TBC-1, which corresponds 
to Tw — To = Tw — Tæ, the second case was TBC-2 corresponds to 
Tw1 — Too = 2(Tw — Too), the third case was TBC-1/2 corresponds to 
Tw1 — Too = (Tw — Too)/2, and the fourth case was TBC-A corresponds 
to adiabatic upstream wall and step (where Ty is the downstream 
wall temperature, Ty; is the upstream wall and step temperature, 
and T» is the free stream temperature). It was found that the recir- 
culation region remains attached to the heated wall downstream 
of the step when the upstream wall and the step are heated to 
the same or different uniform temperatures as the downstream 
wall. The detach of the recirculation region from the heated wall 
downstream of the step may occur at high buoyancy levels when 
the upstream wall and step are maintained adiabatic. For higher 
upstream heating levels, it was found that the reattachment length 
decreases as the buoyancy level increases. Furthermore, it was 
noticed that when the upstream wall and the step are heated to 
a higher temperature than the downstream wall, the velocity dis- 
tribution downstream of the step may exhibit an overshoot over its 
free stream value. It was also noticed that the heat transfer from 
the downstream wall to the fluid decreases as the upstream wall 
and the step are heated. 

Iwai et al. [38] presented a numerical study to investigate the 
effects of buoyancy-assisting laminar mixed convective flow and 
thermal fields through a 3D vertical duct having backward-facing 
step. The downstream wall was maintained to a uniform heat flux, 
while the rest ofthe walls were assumed to be adiabatic. The simu- 
lations were carried out for buoyancy parameter range between 
0<£&<0.12, other parameters such as Reynolds number, aspect 
ratio, expansion ratio, were maintained to be fixed at 125, 16, 2, 
respectively. Iwai and his coworkers have noticed that 2D are effec- 
tive in predicting such type of phenomena for buoyancy parameter 
& « 0.06, for values greater than this, a 3D computational study 
becomes necessary. The reattachment point and the peak Nusselt 
number move closer to the step as the buoyancy force increases. 
Moreover, the secondary recirculation region which developed at 
the corner between the step and the stepped downstream wall 
increases in size as the buoyancy force increases. It was found that 
the maximum Nusselt numbers are located symmetrically near the 


side walls far from the center plane of the downstream wall. In 
addition, the maximum Nusselt number was noticed to be down- 
stream of the reattachment point and the distance between them 
increases with increasing buoyancy force. 

The research to study effects of separation and reattachment 
phenomenon does not concern on laminar flow. Abu-Mulaweh 
et al. [39] examined the effect of turbulent mixed convection over 
a 2D vertical backward-facing a step. The step geometry consists 
of an adiabatic backward-facing step, an upstream wall and a 
downstream wall heated at constant and uniform temperatures. 
A free stream air was used as a medium in the velocities range of 
0 < u% x 0.41 m/s fora step height of 22 mm, and a temperature dif- 
ference (AT- 30?C) between the heated walls and the free stream. 
The measurements for different free stream velocities (us; =0, 0.23, 
0.41 m/s) were carried out at one streamwise location upstream of 
the step (x » —5 cm) and four different streamwise locations down- 
stream of the step (x 3.5, 6, 10, 25 cm). It was indicated that the 
streamwise velocity increases as the free stream increases and 
vice versa for the temperature distribution for both upstream and 
downstream locations. The reattachment length was noticed to be 
increased with increasing the free stream velocity. The magnitudes 
for the intensities of streamwise and transverse velocity fluctua- 
tions were found to be larger for natural convection than for mixed 
convection on the upstream location. While, the intensities tem- 
perature fluctuations decrease rapidly for mixed convective flow. 
On the other hand, the magnitudes of the intensities of both veloc- 
ity and temperature fluctuations are smaller for mixed convection 
than natural convection at the four downstream locations. It was 
found that the reattachment length increases while the Nusselt 
number from the downstream heated wall decreases as the free 
stream increases, and as a result the location of the maximum heat 
transfer rate moves away from the step as the free stream velocity 
increases. 

Abu-Mulaweh et al. [40] studied experimentally the effect of 
the backward-facing step heights on the turbulent mixed convec- 
tion flow along a 3D vertical flat plate. The experimental setup, step 
geometry, boundary conditions, and measurement equipments 
were identical to the previous work done by the same authors in 
2001 [39]. The experiments were carried out for step heights of 
0, 11, and 22mm, while the free stream air velocity, us, and the 
temperature difference between the heated walls and free stream, 
AT, were fixed at 0.41 m/s and 30°C, respectively. Measurements 
of the flow field were carried out at a location upstream of the 
step (x=—5 cm) at the midplane (z- 0) of the plate's width for pure 
forced, pure natural, and mixed convection flows. It was found that 
there was a velocity overshoots in the mixed convective case due 
to the buoyancy assisting near to the heated wall upstream of the 
step for the dimensionless streamwise velocity distributions. It was 
Observed that increasing the step height enhances the turbulence 
intensity, which causes the flow to become turbulent downstream 
of the step. This is because the existence of a backward-facing step 
triggers transition from the re-laminarized flow upstream of the 
step to turbulent flow downstream of the step. Furthermore, it was 
found that as the step height increases, the turbulent intensity of 
both streamwise and transverse velocity fluctuations, the intensity 
of the of temperature fluctuations, and heat transfer rate down- 
stream of the step increase. The length of the recirculation region 
increases with increasing the step height, which results the loca- 
tion of the maximum heat transfer rate moves away from the step 
as the step height increases. 


4. Fundamentals of nanofluids 


Heat transfer is commonly encountered in many engineering 
systems for heating and cooling. The high performance cooling is 
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essential in many industrial systems such as nuclear plants, HVAC, 
transportation, energy production, supply to electronics, and paper 
production. It is also important in consumer products such as com- 
puters, car engines, power electronics, and high-powered lasers or 
X-rays. Traditional fluids are used as coolant and their heat transfer 
increase with thermal conductivity increment. The development 
in technologies to perform desired performance, reliable, and min- 
imize size of products reveal to increase in heat loads which in 
some cases it exceeds 25 kW, therefore liquids and gases have lim- 
ited thermal performance due to their low thermal conductivity. 
Many researchers tried to increase the thermal conductivity of 
the conventional fluids by suspending millimeter- or micrometer- 
sized solid particles in liquids since the thermal conductivity of 
solids are higher than liquids. Maxwell [41] made the first effort 
using millimeter-sized particles suspended in water, but due to 
the large size and high density the small particles settled out of the 
suspension. 

Many trails have been done to disperse solid particles in liq- 
uids since Maxwell's time and all of them observed settlement due 
to the large size and density of the particles. In 1993, Choi con- 
ceived the novel concept of dispersing nanometer-sized particles 
into base fluids and he proposed the name of Nanofluids in 1995 
[42]. Nanofluids are dilute suspensions of nanoparticles smaller 
than 100 nm, which belong to a new type of functional composite 
materials. These solid-liquid composites are very stable and show 
higher thermal conductivity and higher convective heat transfer 
performance than traditional liquids. Nanofluids have small size 
which allows them to have less particle sedimentation, erosion, 
passage clogging, and pressure drop than microfluids or fluids with 
larger scaled suspensions. 

The interest of studying nanofluids have been increased since 
Choi's concept is presented, and many researches were conducted 
in different fields and facilities to understand the structures, func- 
tions, mechanisms, applications, environmental impacts, and other 
characteristics of nanofluids. There were some researchers made 
review papers such as Yu et al. [43], Wang and Mujumdar [44], Li 
et al. [45], Kakac and Pramuanjaroenkij [46], Wen et al. [47], Das 
et al. [48], Wong and Castillo [49], and Nsofor [50]. They declared 
that nanofluids will produce a new generation of cooling technol- 
ogy, and further researches is needed to improve the stability to 
avoid cluster phenomena in nanofluids. 


4.1. Production of nanoparticles 


The invention and improvement of modern technologies per- 
mits to fabricate materials with nanometer scale. Many types 
of materials have been used to make nanoparticles that used in 
nanofluids. Nanoparticles are synthesized by two different pro- 
cesses which are physical process and chemical process. In physical 
process, mechanical grinding methods using grinding balls, and 
the inert-gas-condensation technique were utilized. However, 
in chemical process many methods are used such as chemical 
precipitation, chemical vapor deposition, micro-emulsions, spray 
pyrolysis, and thermal spraying. Carbon nano-tubes (CNTs) and 
fullerenes are commonly formed in ordinary flames produced by 
burning methane, ethylene, and benzene, and highly irregular in 
size and quality in nature. These types of nanoscaled particles are 
formed by different processes such as arc discharge, laser ablation, 
chemical vapour deposition, and ball milling. 


4.2. Nanoparticle material types 


Certain materials have been used to produce nanoparticles and 
disperse them in host liquids such as oxide ceramics (Al,03, CuO); 
nitride ceramics (AIN, SiN); carbide ceramics (SiC, TiC); metals 
(Ag, Al, Au, Cu, Fe); semiconductors (SiO2, TiO); single, double, 


or multi wall carbon nanotubes (SWCNT, DWCNT, MWCNT); and 
composite materials such as nanoparticle core-polymer shell com- 
posites. Moreover, new materials and structures are desired for use 
in nanofluids where the solid-liquid interface is doped with various 
molecules. 


4.3. Host liquid types 


There are many types of metallic and nonmetallic liquids have 
been used as host fluid to produce nanofluids. Metallic liquids 
such as sodium (644 K) have higher thermal conductivity than non- 
metallic liquids such as water and organic liquids (ethylene glycol, 
oil). Water is used as the best traditional liquid due to its high 
thermal conductivity, abundance, low cost, and friendly to the envi- 
ronment. 


5. Preparation of nanofluids 


The powder form nanoparticles which disperse in host liq- 
uids are called nanofluids. Nanofluids can be produced by two 
techniques; the two-step (double-step) method, and one-step 
(single-step) method. These methods have been utilized using dif- 
ferent types of chemical and physical techniques to make sure that 
the solid-liquid mixture is stable to avoid agglomeration, additional 
flow resistance, possible erosion and clogging, poor thermal con- 
ductivity, and poor heat transfer. The two-step method is done 
by producing the nanoparticles powder initially as introduced in 
the previous section, and then disperses them into a host liquid. 
However, in one-step method the nanoparticles are simultaneously 
made and directly dispersed into the base fluid [51]. It is noticed in 
the literature that nanofluids with oxide nanoparticles and carbon 
nanotubes are produced well by the two-step method, while it is 
not suitable for nanofluids with metallic nanoparticles. 


5.1. Techniques of nanofluids production 


5.1.1. Two-step technique 

It is introduced before that the two-step method is started 
by nanoparticles production separately and afterward dispersed 
into a base fluid. The advantage of this method is that it is easily 
and economically produced. On the other hand, the disadvantage 
of this method is the quick agglomeration of individual particles 
before the achievement of complete dispersion due to Vander 
Waals attractive forces between nanoparticles. This agglomera- 
tion is a big obstruction to achieve high heat transfer performance 
due to the quick settling of particles out of the base fluids 
and it becomes worse as the volume concentration increases. 
The agglomeration is not only a problem in nanofluid technol- 
ogy, but it is also a critical issue in all nanopowder technology 
especially during drying, storage, and nanoparticle transportation 
stage. 

Many different methods have been used to avoid agglomeration 
using two-step process to initiate the move towards commercial- 
ization by facilitating the mass production of nanofluids. Two-step 
methods such as stirrer, ultrasonic bath, ultrasonic disrupter, and 
high pressure homogenizer are the popular methods of two-step 
technique. Eastman et al. [52], Lee at al. [53], and Wang et al. [54] 
used ultrasonic equipment to disperse Al203 in the base fluid with 
less agglomeration. Moreover, Hong et al. [55,56] utilized ultrasonic 
equipment produced Fe nanofluids by mixing Fe nanocrystalline 
powder in ethylene glycol. Ultrasonic equipment was also used 
to produce TiO;-water nanofluid as introduced in Murshed et al. 
[57]. Metallic nanofluids were produced by using different two-step 
methods, but they were not successful due to the high agglom- 
eration. There are some other techniques that are used to attain 
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the stability of the nanoparticles in the base fluids using the two- 
step methods such as controlling pH, or adding surfactants and 
dispersant. These techniques are explained in detail in the next 
section. 


5.1.2. One-step technique 

One-step technique is the process of producing nanofluids 
where the nanoparticles are simultaneously made and directly dis- 
persed into the base fluid. This method is preferable to produce 
nanofluids containing high thermal conductivity metals to avoid 
erosion and oxidation of particles. The advantage of this process 
lies on minimizing of nanoparticles agglomeration. This behavior 
increases the stability of the suspensions and uniform dispersion 
in the host liquids. The disadvantage of one-step technique is the 
limit of quantity ofthe production due to the slow ofthe production 
process, low concentration of nanoparticles, and the high cost. 

Different methods have been used to reduce the time and cost 
using evaporation, physical or chemical one-step methods. The 
first approach of one-step evaporation method is direct evapora- 
tion method which called Vacuum Evaporation onto a Running Oil 
Substrate (VEROS) technique developed by Akoh et al. [58]. This 
method was created to produce nanoparticles, but it was found 
that it is difficult to separate the particles from fluids to produce 
dry nanoparticles. Later, Wagener et al. [59] proposed a modified 
VEROS technique using high pressure magnetron sputtering to pro- 
duce metallic nanofluids such as Ag and Fe. Furthermore, Eastman 
et al. [52] modified VEROS technique by using direct condensed 
Cu vapour to form nanoparticles and in contact with flowing with 
low-vapor pressure ethyl glycol (EG). The well dispersed of Cu 
nanoparticles in ethyl glycol enhanced the thermal conductivity of 
the base fluid by up to 40% at the particle volume concentration of 
0.3 vol.%. For one-step physical method; a submerged arc nanopar- 
ticles synthesis is used in [60-62] to produce nanofluids of Cu, CuO, 
and TiO2. In this method, nanoparticles are produced by heating 
the solid material from an electrode in the form of arc sparkling 
and condensed into liquid in a vacuum chamber to form nanofluid. 
Finally, the one-step chemical method is reported by Zhu et al. [63]. 
They produced Cu nanoparticles of less than 20 nm and dispersed 
in (EG) by the reduction of a copper salt by sodium hypophosphite 
under microwave irradiation. A protective polymer and stabilizer 
called polyvinylpyrrolidone was added so that it inhibited agglom- 
eration. Wei et al. [64] developed a chemical solution method (CSM) 
to synthesize Cu5O nanofluids by suspensions of cuprous-oxide 
(Cu50) in water. Their experimental studies showed a substantial 
thermal conductivity enhancement of up to 24% with their synthe- 
sized nanofluids. Among all the three methods; chemical process 
is the fastest and can be produced in large quantities but it is still 
not in commercial quantity as in two-step process. 


5.2. Stability 


The stability of the suspensions in nanofluids is the most criti- 
cal issue for enhancing heat transfer. This is because nanoparticles 
tend to aggregate with the time elapsed for its high surface-activity, 
which results not only the settlement and clogging of microchan- 
nels but also the decreasing of thermal conductivity of nanofluids. 
Nanofluids are considered to be stable when the concentration 
or particle size of supernatant particles keeps constant. There are 
some methods have been used to determine the stability of nanoflu- 
ids such as sedimentation method [65,66], Zeta potential method 
and absorbency method [67]. Sedimentation method is the most 
simple and reliable method. It is basically obtained by special appa- 
ratus which functioning to measure the variation of concentration 
or particle size of supernatant particle with sediment time. More- 
over, another method using an optical sedimentation photograph 
to monitor the stability of nanofluids was used but its disadvantage 


is the long time of the observation compared to other methods. On 
the other hand, Zeta potential method has a limitation ofthe viscos- 
ity and concentration of nanofluids. A UV-vis spectrophotometer 
was used to detectthe stability of nanofluids. The variation of super- 
natant particle concentration of nanofluids with sediment time 
can be obtained by the measurement of absorption of nanofluids 
because there is a linear relation between the supernatant nanopar- 
ticle concentration and the absorbance of suspended particles. It is 
declared that the most important factors affecting the stability of 
suspensions were the nanoparticles concentration, additives, vis- 
cosity of base liquid and pH value. The variety, diameter, density of 
nanoparticle and ultrasonic vibration also influence the stability of 
nanofluids. In addition, the particle morphology, chemical structure 
of particles and base fluids affect the stability as well. It is observed 
that the stability of nanofluids do not exceed more than few weeks 
or months which need to have new techniques to improve the time 
of nanofluids stability. 


6. Experimental investigations on nanofluids 
6.1. Thermal conductivity 


Itis reported that many experimental studies have been done on 
thermal conductivity of nanofluids because it is a primary assess- 
ment of the heat transfer performance of nanofluids. The physical 
mechanism for thermal conductivity enhancement of nanofluids is 
not well understood up to date. Maxwell [68] was one of the first 
investigators who suspended particles thermal conductivity. He 
used a very dilute suspension of spherical particles by ignoring the 
interactions between the particles, and made an equation that can 
be applied only to mixtures with low particle volume concentra- 
tions. The gap of particles interactions attracts many researchers to 
extend or improve the Maxwell equation with various factors that 
affect thermal conductivity such as particle shape [69-77], parti- 
cle dispersion and distribution [78], volume concentration [79-82], 
particle-shell structure [69,83-88], and contact resistance [89,90]. 
Even though these studies tried to improve equations to predict 
thermal conductivity but the comparison between the experimen- 
tal data and theoretical predictions for nanofluids is generally not 
satisfactory. New techniques of mechanism have been used to 
improve the prediction of nanoscale such as nanoparticle-matrix 
interfacial layer [91-94], nanoparticle Brownian motion [95-99], 
and nanoparticle cluster/aggregate [100,101]. 


6.1.1. Thermal conductivity measurement methods 

There are many methods used to measure nanofluids thermal 
conductivity such as transient hot-wire (THW), temperature oscil- 
lation (TO) and steady-state (SS) methods [55-57,61,102-105]. 
Transient hot-wire method is well established as the most accu- 
rate, reliable and robust measurement technique for the thermal 
conductivity of nanofluids, and so it is utilized in most of thermal 
conductivity measurements. The transient hot-wire is functioned 
as both heater and thermometer and many modifications have been 
proposed such as coating the hot wire. In oscillation method; the 
electrical components are removed from the test sample and per- 
form as purely thermal as reported by [106-108], so that the ion 
movement is not affecting the measurements. The steady-state 
method was used by [54,109] to measure the effective thermal 
conductivity of Al;05 and CuO nanofluids with an improvement 
of 12% in thermal conductivity enhancement at 3% volume fraction 
of nanoparticles. 


6.1.2. Effects of thermal conductivity parameters 

Many experimental studies have been done by different 
researchers for examining and improving the thermal conduc- 
tivity enhancement of certain types of nanofluids using different 
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methods and parameters. They calculated the thermal conductivity 
enhancement by using the ratio of the thermal conductivity of the 
nanofluid over its base fluid thermal conductivity. It is noticed that 
most of the studies did not concerned or noticed the particle size 
distribution. In addition, the effect of the particles agglomeration 
in the base fluid is also unknown. 

It is declared that the thermal conductivity of nanofluid is sen- 
sitive to certain parameters which are responsible for constant 
nanofluids such as particle volume concentration, particle size, 
particle material, particle shape, base fluid material, temperature, 
additives and surfactants, and acidity. These parameters will be 
introduced in the following section separately base on the previous 
experimental investigations and findings. 


6.1.2.1. Effects of particle volume concentration. One of the most 
important parameters is volume concentration or in other words 
volume fraction. It is the ratio of the volume of nanoparticles to 
the base fluid volume. As declared previously that solids have 
higher thermal conductivity than liquids, so as the percentage 
of nanoparticles increases, the thermal conductivity increases as 
well. Unfortunately, this is not working for all cases; it is declared 
that many colloidal and biological suspensions show strong non- 
Newtonian behavior but nanofluids with volume concentration not 
more than 5% exhibit Newtonian behavior. At higher particle vol- 
ume concentrations, the increase in enhancement is expected to 
diminish or even reverse, and nanofluid viscosity would increase. 
Three methods have been used by Wang et al. [54] to measure the 
viscosity of nanofluids but they did not observe any non-Newtonian 
effects. They used Al205-water nanofluid and found a 30% increase 
in viscosity when compared with pure water at 3% volume con- 
centration of the particles. However, Pak and Cho [110] measured 
the viscosity to be much higher. There are many researches done 
on particle volume concentration [42,54,102,106,109,111-126] for 
different particle size, particle shape, particle material, base fluid 
material, additive, and acidity level. They reported that at low vol- 
ume concentration (<5%) the thermal conductivity ratios increase 
linearly with volume concentration. Furthermore, it is noticed that 
metallic nanofluids and either single- or multi-wall nanotubes 
enhance thermal conductivity with low volume concentration 
compared to oxides and nonmetallic nanofluids. 


6.1.2.2. Effects of particle size. The size of particle plays a significant 
role in thermal conductivity and heat transfer enhancement in base 
fluids. As introduced before that different suspension sizes of mil- 
limeter and micrometer particles in host liquid have been used to 
increase the thermal conductivity, but the particles agglomerate 
quickly and settled out of the liquid. Nanoparticles are formed to 
decrease the size of the particles so that the time of sedimenta- 
tion decreases to reach in some cases to more than few weeks or 
months. The size of nanoparticles are in the range of 1-100 nm and 
as the size decreases the thermal conductivity increases due to the 
increase of the relative surface area [127]. The much larger relative 
surface area of nanoparticles, should not only significantly improve 
thermal conductivity capabilities, but also should increase the sta- 
bility ofthe suspensions. In addition, it is reported that as the size of 
nanoparticle reduced, the Brownian motion will be induced [128]. 
Furthermore, there are some researches indicate that the thermal 
conductivity enhanced with nanoparticles increment to ~60nm 
and vice versa for greater values [53,114]. The size of nanopar- 
ticles in the host liquid is too important to the modern science 
and technology. It can be applied in small size systems such as 
micro-channels and nanosystems where the clog and erosion is 
prohibited. 


6.1.2.3. Effects of particle material. There are many types of materi- 
als have been utilized to produce nanofluids such as oxide ceramics 


(Al205, CuO); nitride ceramics (AIN, SiN); carbide ceramics (SiC, 
TiC); metals (Ag, Al, Au, Cu, Fe); semiconductors (SiO2, TiO;); sin- 
gle, double, or multi wall carbon nanotubes (SWCNT, DWCNT, 
MWCNT); and composite materials. It was found that metallic 
and allotrops of carbon nanoparticles (such as diamond) increase 
the thermal conductivity of a host fluid more than other types of 
nanoparticles having the same volume concentration. 


6.1.2.4. Effects of particle shape. Thermal conductivity enhance- 
ment in nanofluids is affected by the geometrical shape of the 
particles. Different shapes such as spherical and cylindrical (fiber, 
rod, and tube) nanoparticles have been used in different host liquids 
to identify the enhancement of nanofluids thermal conductivity. It 
is reported that cylindrical nanoparticles are better than spherical 
shape in enhancing the thermal conductivity of nanofluids [57,129]. 
This is due to the elongated particles and the large aspect ratio of 
the surface area to volume that conducts heat through the fluid. 


6.1.2.5. Effects of base fluid material. Many host liquids have been 
used to produce nanofluids such as aqueous and organic liquids 
(ethylene glycol, and oils) to enhance thermal conductivity. It is 
found that the thermal conductivity enhancement increases as the 
thermal conductivity of the base fluid decreases [114,115,130]. 


6.1.2.6. Effects of temperature. The nature of the nanofluid ther- 
mal conductivity is sensitive to temperature. Many studies have 
been done for the effect of temperature on thermal conductiv- 
ity. They reported that the thermal conductivity enhanced as the 
temperature increases due to the motion of the nanoparticles 
[106,109,127,131]. In opposite, Masuda et al. [111] found that the 
thermal conductivity enhancement decreases as the temperature 
increases. 


6.1.2.7. Effects of additives and surfactants. Itis declared that certain 
types of additive have been used in nanofluids to keep nanoparti- 
cles in suspension form and prevent agglomeration for improving 
the thermal conductivity enhancement. The additives types and 
concentration is added depend on the type of both nanoparticles 
and base fluids. Eastman et al. [113] used thloglycolic acid in Cu- 
ethyl glycol, Assael et al. [132] used CTAB and Nanosperse AQ for 
MWCNT-water, Ding et al. [122] gum Arabic for MWCNT-water, 
and other types of acidic fluids such as Oleic acids, etc. 


6.1.2.8. Effects of Acidity (PH). Limited studies have been published 
on the effect of fluid acidity on the thermal conductivity enhance- 
ment of nanofluids [117,126]. But the general trend is that acidity 
increases the thermal conductivity enhancement. 


6.1.3. Mechanisms of thermal transport in nanoscale 

6.1.3.1. Brownian motion. Brownian motion is one of the param- 
eters that affect the flow and heat transfer of nanofluids. It is 
basically the random dynamic mode of particles in a liquid where 
the particles collide between each other. It was determined that 
the thermal conductivity increases with an increase in tempera- 
ture, but also it was shown that nanofluids compose of smaller 
particles experienced greater enhancement than with larger par- 
ticles. Since temperature represents the overall kinetic energy of 
the particles, an increase in temperature will cause increase in 
the particles motion. It is easier for smaller particles to move; 
therefore, smaller particles will display a higher level of Brown- 
ian motion than larger particles. It is outlined that the movement 
of nanoparticles due to Brownian motion was too slow in trans- 
porting heat through a fluid and it can be ignored [133-136]. To 
travel from one point to another, a particle moves a large dis- 
tance over many different paths in order to reach a destination 
that may be a short distance from the starting point. Therefore, 
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the random motion of particles, no matter how agitated or ener- 
getic they may be not a key factor in the improvement of heat 
transfer. On the other hand, there are some studies developed a 
dynamic model that takes into account convection heat transfer 
induced by Brownian nanoparticles. They reported that Brownian 
motion of nanoparticles at the molecular and nanoscale levels is a 
main mechanism controlling the thermal conductivity of nanoflu- 
ids [96,98,137]. 


6.1.3.2. Liquid layering on the nanoparticle-liquid interface. There 
are some studies declared that there is a liquid layering on the 
nanoparticles, which improves the enhancement of nanofluid heat 
transfer properties. The thickness and thermal conductivity of the 
nanolayer are not known yet, but the liquid molecules close to a 
solid surface have been proven by Yu et al. [138] to form layers. 
Ren et al. [94] indicated that the thermal conductivity increases 
with liquid layer increment. They assumed that the thermal con- 
ductivity of the layer would be somewhere between the thermal 
conductivity of the bulk fluid and the nanoparticle. On the other 
hand, Keblinski [135] assumed that the thermal conductivity of the 
interfacial liquid is the same as that of the solid and showed that 
the resultant large effective volume of the particle-layered-liquid 
structure would enhance the thermal conductivity. Ren et al. [94] 
also found that as the nanoparticles increased in size, the effects of 
the liquid layering became weaker. Yu and Choi [92] determined 
that the solid-liquid interfacial layers in the nanofluid plays an 
important role in the enhanced thermal conductivity especially 
when the particle diameter is less than 10 nm. 


6.1.3.3. Nature of heat transport in nanoparticles. Heatis transferred 
via conduction implies that heat is transferred internally, by vibra- 
tions ofatoms and molecules. The vibrations ofthe atoms which are 
connected together as give rise to the vibrations of the whole crys- 
tal, that is, to lattice vibrations. The energy of the whole vibrating 
system is quantized, and the quantum of thermal energy absorbed 
or emitted by an atom is called a phonon. Therefore, a phonon is 
essentially a quantized mode of vibration occurring in rigid crys- 
tal lattice and plays a major role in material's thermal conduction. 
Keblinski [135] demonstrated that the phonon mean free bath is 
much shorter in the liquid than in the particle and would be effec- 
tive as the distance between the particles and the thickness of the 
layered liquid are small. A greater phonon density exists in hot 
regions of a crystal than in cooler regions. Therefore, heat conduc- 
tion is essential due to the diffusion of phonons, and because the 
temperature gradient changes from hot to cold regions. Electrons 
can also carry heat; thus, metals have many free electrons which 
move around randomly, heat can be transferred from one part ofthe 
metal to another quite effectively. That is why metals are generally 
very good conductors of heat. 


6.1.3.4. Effects of nanoparticles clustering. Clustering of particles is 
one of the major problems that affects the enhancement of ther- 
mal conductivity by creating path of lower thermal resistance. 
It is noticed that as the particles agglomerate increases, the size 
increases and becomes denser where they clustered down, and 
end with low thermal conductivity enhancement. To avoid or min- 
imize the particles clustering, nanofluids should consist of low 
volume concentration and small particles size. The effective vol- 
ume of a cluster, which is the volume from which other clusters 
are excluded, can be much larger than the physical volume of the 
particles [135,136]. The effective volumetric concentration of the 
particles, inside a cluster, is larger than the volume of the solid 
phase, since within such cluster heat can move very rapidly. 


6.1.3.5. Thermophoresis. Thermophoresis refers to the motion of 
colloidal particles in response to a temperature gradient which 


can be clarified by applying the kinetic theory. The high energy 
molecules in the hot region of the liquid impinge on the particles 
with greater momentum than do molecules coming from the cold 
region, thus leading to the migration ofthe particles in the direction 
opposite to the temperature gradient. This mechanism has been 
found to be negligible, primarily because of the very small particle 
migration velocity [91,136]. 


6.1.3.6. Reduction in thermal boundary layer thickness. A few num- 
bers of researches have been done for studying the thermal 
boundary layer such as in Ding et al. [139]. They mentioned 
that a reduction in the thermal boundary layer thickness may 
be a mechanism that causes heat transfer enhancements in 
nanofluids. It is still needed to have more research on this mech- 
anism to understand its effects in nanofluids convective heat 
transfer. 


6.2. Viscosity 


Nanofluids are expected to be used under flow conditions and 
the flow of suspension which are sometimes drastically different 
from that of most common heat transfer fluids that have Newtonian 
characteristics, so itis essential to have the rheological properties of 
nanofluid to use practically. Rheological properties can provide the 
knowledge on the microstructure under both static and dynamic 
conditions, which are important to understand the mechanism of 
heat transfer enhancement using nanofluids. Fluids have frictional 
forces between the molecules and, therefore, they display a certain 
flow resistance which can be measured as viscosity. The fluids are 
either Newtonian or non-Newtonian, depends on their viscosity. It 
has been demonstrated by Newton that the shear force acting on 
a liquid is proportional to the resulting flow velocity. If the viscos- 
ity of the fluid remains constant with an increase in shear rate, it is 
called Newtonian fluid, but if the viscosity changes with an increase 
in shear rate, it is called non-Newtonian fluid. The shear viscosity 
of samples displaying shear-thickening behavior is dependent on 
the degree of shear load. However, the viscosity increases with an 
increase in shear stress. It is declared that the viscosity of nanoflu- 
ids is mainly dependent on volume concentration, nanoparticle 
size, and temperature [140-148]. Nanofluids with high concen- 
tration and large size nanoparticles have high particle interaction 
and may result in particles agglomeration and thus increase the 
flow resistance and decrease the heat transfer enhancement. Par- 
ticle shape plays an important role since during the shear process, 
particles rotate as they move. Cube-shape particles take up more 
volume when rotating than spherical particles and hence less free 
volume is available for the liquid between the particles. In addi- 
tion, it was reported that the viscosity decreases with temperature 
increment. 


6.3. Heat capacity 


The heat capacity of a substance is the amount of heat required 
to change its temperature by one Kelvin, and has units of joule per 
Kelvin (J/K) in the SI system. The specific heat capacity is essen- 
tially the heat capacity per unit mass, usually per gram of material 
(J/g K). The specific heat capacities of nanofluids are different from 
that of base fluid and increase with the size and volume concen- 
tration of nanoparticles decrement. The high specific interfacial 
area of nanoparticle can absorb liquid molecules to its surface and 
form liquid layers, which will reversely constrain nanoparticle and 
turns its free boundary surface atoms to be non-free interior atoms. 
This effect will also increase the specific heat capacity of nanofluid 
[149-153]. 
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7. Applications of nanofluids 


Nanofluids could be used as a super coolant in many industrial 
and consumer products such as nuclear reactors, car engines, radi- 
ators, computers, X-rays, and many other applications. It is called 
super-coolant because it can absorb heat more than any traditional 
fluids, so it can minimize the size of system and increase its per- 
formance. Most of the researches were concerned on studying the 
properties, parameters, and heat transfer enhancement, but not 
much research done for commercial systems. Some researchers 
tried to apply nanofluids in real systems and observe their flow 
and heat transfer. Tzeng et al. [154] studied the effect of nanofluids 
when used as engine coolants. CuO (4.4 wt%) and Al203 (4.4 wt%) 
nanoparticles and antifoam were individually mixed with auto- 
matic transmission oil, and experiment them by a real drive of 
four-wheel-drive (4WD) transmission system. The experimental 
results showed that, the antifoam-oil provided the highest tem- 
perature distribution in rotary blade coupling and, accordingly, the 
worst heat transfer effect, and CuO-oil provided the lowest tem- 
perature distribution both at high and low rotating speed and, 
accordingly, the best heat transfer effect. Chein and Huang [155] 
used Cu nanofluid with various concentrations to study silicon 
microchannel heat sink performance. They found that nanofluid 
cools better than pure water and observed no extra pressure drop 
because of small particle size and low particle volume fraction. 
Koo and Kleinstreuer [156] used two types of nanofluids (i.e., CuO 
nanospheres at low volume concentrations in water and in ethy- 
lene glycol) for the conjugate heat transfer problem for microheat 
sinks. The effect of Brownian motion on the effective fluid viscos- 
ity was considered and found to be less significant than that on 
the effective thermal conductivity. Many other industrial and con- 
sumer companies have their own researches using nanofluids such 
as medical equipments, computer manufacturing, and gas and oil 
[157]. 


8. Safety 


Up to date the behavior, parameters and mechanisms of 
nanofluids are not clearly understood, and because of these vari- 
ables scientists still cannot accurately predict how nano-materials 
will affect living organisms. Now, it is just clear that the biolog- 
ical activity and biokinetics of nanoparticles are different from 
larger particles, and that they depend on many parameters. These 
parameters can modify cellular uptake, protein binding, transloca- 
tion from portal of entry to target site, and the possibility of tissue 
injury. Several studies have demonstrated that nanoparticles are 
more toxic than micro-particles of the same material. 


9. Environmental effects 


Nanofluids are observed to be an environmental friend due to its 
high heat transfer enhancement. The usage of nanofluids decreases 
the consumption of fuel in engines, and electricity which are the 
highest pollutants that threat our earth and increase its temper- 
ature. It is noticed that if nanofluids improve chiller efficiency by 
1%; it will save 320 billion kwh of electricity or equivalent to 5.5 
million barrels of oil per year just in USA alone. 


10. Convective heat transfer of nanofluids 


Convective heat transfer is the phenomena of the macroscopic 
motion of the fluid relative to the surface due to the heat trans- 
fer between them. The surface could be a solid wall or an interface 
with another liquid. Convective heat transfer is divided into two 
major types, which are natural convective heat transfer where fluid 


motion is induced by buoyancy, and forced convective heat transfer 
where fluid is forced to flow through a boundary region. Convective 
heat transfer of a medium depends on its thermophysical proper- 
ties such as, thermal conductivity, specific heat capacity, viscosity, 
density, and thermal expansion coefficient. 


10.1. Natural convection 


The natural convection studies using nanofluids are very few. 
Even though the research is limited, however the available results 
are diverged. Putra et al. [158], and Wen and Ding [120] made 
an experimental investigation for natural convection using Al203 
and TiO; nanofluids, respectively. They found that the presence of 
nanoparticles in water systematically decreased the natural con- 
vective heat transfer coefficient. Furthermore, they declared that 
the decrease was due to the effects of particle/fluid slip and sedi- 
mentation of nanoparticles. In contrast, Khanafer et al. [159] used a 
numerical technique and predicted that nanofluids enhanced nat- 
ural convective heat transfer. It is needed to do more research on 
natural convection to get a clear view and well understanding for 
such phenomena. 


10.2. Forced convection 


Forced convection studies using nanofluids are also limited and 
they have diverged results from each other. Majority of forced con- 
vective heat transfer using nanofluids found that there is convective 
heat enhancement such as Li and Xuan [160], Xuan and Li [161], 
Lee and Choi [162], Jang and Choi [163], and Xuan and Roetzel 
[164]. There are some studies reported that heat transfer can be 
enhanced under certain conditions (such as low volume concen- 
tration for decreasing the viscosity) by Pak and Cho [165], Lee and 
Mudawar [166], and Chein and Chuang [167]. A very few studies 
show a decrease in the convective heat transfer coefficient when 
nanoparticles are added to the base liquids such as Yang et al. [168]. 


10.3. Convection of nanofluids in BFS 


It is noticed in the literature that there are few experimen- 
tal and numerical studies have been done to investigate the flow 
and heat transfer of nanofluids in different geometries with dif- 
ferent regimes. Engineering geometries such as tubes and annuli 
were the first geometries that have been used for investigating 
nanofluids enhancement. Most of the studies show that nanoflu- 
ids have higher thermal transport than any ordinary fluids and 
their enhancement increase or decrease depends on the change 
of the parameters and characteristics which introduced in the pre- 
vious sections. Other studies have been done to apply nanofluids 
for microchannels, different geometrical shapes, and systems. The 
results show an improvement in thermal transport enhancement. 

To the best knowledge of the authors, there is only one paper 
investigated the flow and heat transfer influence of nanofluids 
through a duct with backward facing step. Abu-Nada [169] stud- 
ied the flow and heat transfer influences over a backward facing 
step (BFS) using nanofluids numerically. Different volume fractions 
and types of nanoparticles were used in the base fluid (water) and 
flowed through BFS channel with an expansion ratio ER=2. The 
Reynolds number and nanoparticles volume fraction used were 
in the range of 200 < Re x 600 and 0<@<0.2, respectively. The 
Prandtl number of the base fluid was kept constant at 6.2 for the 
five types of nanoparticles which are Cu, Ag, Al203, CuO, and TiO». It 
was noticed that at the top wall and for Re = 200, an enhancement in 
the Nusselt number was registered by increasing the volume frac- 
tion of the nanoparticles. Moreover, there was an enhancement 
in Nusselt number at the top wall for Re» 300 except in the sec- 
ondary recirculation zone where insignificant enhancement was 
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Table 2 


Summary of the parameters that affect the flow over backward-facing step. 


As the parameter increases 
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Buoyancy-opposing 


Buoyancy-assisting 


Aspect ratio 


Expansion Ratio 


Prandtl number 


Inclination 


Velocity 


Step height 


Decreases 
Decreases 


Increases 


Increases 


Decreases 


Increases 


Decrease, 0° < $ < 180° Increase, 180° < $ < 360° 
Increases, 0° < $ < 180°Decrease, 180° < $ < 360° 


Increases 


Increases 


Local Nusselt number 
Reattachment length 


Increases 


Increases 


R«2.25, increaseR > 2.25, decrease 


Increases 


Increases 


Increases 


registered. It was found that the high Nusselt number inside the 
recirculation zone is mainly depended on the thermophysical prop- 
erties of the nanoparticles and it is independent on of Reynolds 
number. Both the Reynolds number and thermophysical proper- 
ties of the nanoparticles have an influence on the value of Nusselt 
number outside the recirculation region. However, it was noticed 
that outside the recirculation zone, nanoparticles having high ther- 
mal conductivity (such as Ag or Cu) have more enhancements on 
the value of the Nusselt number, and verse visa in the primary and 
secondary recirculation zones. 

Very recently, the same authors [170-172] investigated numer- 
ically the laminar forced and mixed convection assisting flow over 
horizontal and vertical backward-facing steps in a duct using differ- 
ent nanofluids (such as Au, Ag, Al205, Cu, CuO, diamond, SiO2, and 
TiO2). It was found that a primary recirculation region developed 
after the sudden expansion and it starts to become fully developed 
downstream of the reattachment point. The reattachment point is 
found to move downstream far from the step as Reynolds number 
increases. Moreover, they also inferred that low density nanofluids 
have higher absolute velocity compared to high density nanoflu- 
ids. It is found that a recirculation region was developed straight 
behind the backward facing step which was appeared between the 
edge of the step and few millimeters before the corner which con- 
nects the step and the downstream wall. In a few millimeters zone 
between the recirculation region and the downstream wall; a U- 
turn flow was developed opposite to the recirculation flow which 
is mixed with the unrecirculated flow and travels along the chan- 
nel. Two maximum and one minimum peaks in Nusselt number 
were observed along the heated downstream wall. It is inferred 
that Au nanofluid has the highest maximum peak of Nusselt num- 
ber, while diamond nanofluid has the highest minimum peak in the 
recirculation region. Nanofluids with higher Prandtl number have 
higher maximum peak of Nusselt numbers after the separation and 
recirculation flow is vanished. 


11. Conclusions 


In this article, a comprehensive review of previous efforts is 
presented for different convective flow regimes and heat transfer 
through a duct having backward facing step. The effects of sev- 
eral parameters in geometry, boundary conditions, and types of 
fluids were extensively introduced and investigated. In backward- 
facing step, the Nusselt number increases as the parameter such 
as step height, velocity, Prandtl number, and aspect ratio increase, 
and a vice versa for expansion ratio. It was reported that the 
reattachment length increases as step height, velocity, Prandtl 
number, aspect ratio, and expansion ratio less than 2.25 increase, 
and decreases with increasing expansion ratio greater than 2.25. 
These features are summarized in Table 2. In addition, an extensive 
review for preparation, parameters, mechanisms, characteristics, 
convective heat transfer enhancement, applications, safety, and 
environmental impacts for nanofluids was also reported. It is found 
that nanofluids will solve many technological problems and space 
limitation, and will be used commercially such as other ordinary 
coolants in the future. It was observed from the literature review 
that many work have been done using conventional fluids to study 
the heat transfer and fluid flow characteristics in BFS. Thus, more 
work using nanofluids is needed to investigate the heat transfer 
enhancement over backward facing step. 
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